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Abstract. The study of dynamical systems defined on complex 
networks provides a natural framework with which to investigate 
myriad features of neural dynamics, and has been widely under- 
taken. Typically, however, networks employed in theoretical stud- 
ies bear little relation to the spatial embedding or connectivity 
of the neural networks that they attempt to replicate. Here, we 
employ detailed neuroimaging data to define a network whose spa- 
tial embedding represents accurately the folded structure of the 
cortical surface of a rat and investigate the propagation of ac- 
tivity over this network under simple spreading and connectivity 
rules. By comparison with standard network models with the same 
coarse statistics, we show that the cortical geometry infiuences pro- 
foundly the speed propagation of activation through the network. 
Our conclusions are of high relevance to the theoretical modelling 
of epileptic seizure events, and indicate that such studies which 
omit physiological network structure risk simplifying the dynamics 
in a potentially significant way. 



1. Introduction 

The newly emerging discipline of Network Science provides a general 
framework for representing, modelling and predicting the behaviour of 
complex systems belonging to areas as diverse as social science, biol- 
ogy and information technology [1]. Motivated by the observation that 
most real-world networks fail to conform to the homogeneous Poisso- 
nian degree distribution admitted by Erdos-Renyi random graphs [2J, 
improved network models were constructed (most notably the small- 
world model of Watts and Strogatz [3] and the preferential attachment 
model of Barabasi and Albert that were capable of recovering many 
of the interesting features displayed by real-world network data. Initial 
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investigations into complex networks focussed primarily on the char- 
acterisation of networks in terms of a small number of topological pa- 
rameters; however, more recently, interest has shifted towards under- 
standing the influence of network structure on the dynamic processes 
occurring upon them - see, for example, [5l [6] and references therein. 

Such investigations are particularly relevant to biological systems, 
in which a well-defined network structure is frequently a key feature, 
the evolution and topology of which are presumed to affect relevant 
biological processes. A paradigm for such studies is that of neural sys- 
tems, which have been widely studied in this context with the aim of 
providing a more complete understanding of epilepsy and other neural 
conditions [7]. Epilepsy is characterised by recurrent and unpredictable 
instances of "excessive, or synchronous neuronal activity" (seizures) [8J; 
the synchronisation of neuronal activity in networks has therefore re- 
ceived significant attention, and network topology is considered to be a 
dominant factor affecting spreading dynamics [HI [IHl [HI [121 [131 [E], 
independent of the specific model by which activation is transmit- 
ted across the network. A wide range of transmission models has 
been employed in the literature. Representative examples include sim- 
ple cellular automata-type spreading rules [HI [15] and pulse-coupled- 
oscillators [9] (and references therein), of which integrate-and-fire mod- 
els [T3l[T6l|T7] are a special case. We note that in, e.g. [9l[T0l[Tl] (and in 
the present work), network nodes are to be interpreted as 'neural units' 
comprising many neurons, and representing a single cortical column, 
say. Detailed consideration of synaptic signalling models is therefore 
not appropriate. 

A further essential aspect of neural systems is that they are spatially 
embedded - i.e., their nodes and edges are constrained to lie on a fixed 
geometric structure, and it is expected that these additional factors will 
further influence network dynamics. Indeed, though spatial networks 
have received considerable attention of late pS|, the influence of spatial 
embedding upon network dynamics is not well understood. 

A small number of recent studies have begun to employ neuroimag- 
ing data {e.g. via Diffusion Tensor Imaging) to infer large-scale network 
connectivity |19l |20j; however, the majority of studies in macroscale 
epilepsy and seizure modelling typically employs uniform lattices in 
ID or 2D with network connectivity restricted to nearest-neighbours, 
or with additional long-range connections obeying arbitrary rules [211 
|22| [23] . Such an approach fails to exploit a wealth of neuroimaging data 
|24j . which reveals the intricate connectivity structure and surface ge- 
ometry of the brain. Therefore, current approaches still observe simpli- 
fied 2D surfaces [2S1 [2S1 [22] inspired by the original studies of excitable 
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media [28] — an organization that is different from the rounded shape 
of model organisms such as rodents or convoluted brain surfaces such 
as for humans. Note that this limitation not only relates to epilepsy 
modelling but also to studies of spreading depression [291 EO] • 

In this work, we seek to address the deficiency in neural network 
studies discussed above. We employ detailed spatial information ob- 
tained from freely-available neuroimaging data of a rat brain to de- 
fine a physiologically-relevant neural network and, via comparison with 
commonly-employed network architectures, determine the influence of 
the connectivity and complex surface geometry of the brain on seizure 
dynamics. We restrict attention to lattice-like network structures {i.e., 
networks that display significant clustering and long average path- 
lengths), and a basic spreading model in favour of the more complex 
signalling models outlined above, as the aim here is to highlight the 
contribution of spatial network properties to signal transmission. We 
also limit our current study to spreading on the brain surface; the role 
of long-distance white matter fibre tracts between brain regions is not 
part of this work. 

Our investigations highlight clearly that employing cortical surface 
geometry to inform network structure influences dramatically the prop- 
agation of activation through the network. By doing so, we indicate 
that activation dynamics of relevance to epileptic seizure initiation and 
progression (in particular, those highlighting the importance of the site 
of initiation within the network) display distinct differences in idealised 
networks which are commonly employed in the literature. Most strik- 
ingly, we show that our cortical network delays significantly the total 
activation of the network when compared to idealised networks with 
the same coarse statistics: in the parameter regimes that we study, 
the time to activation is increased by a delay factor lying in the range 
1.45-1.88. In this way, we highlight clearly the importance of realis- 
tic network structure on activation dynamics and indicate that such 
considerations should be included in theoretical models which aim to 
provide a more complete description of (e.g.) epileptic seizure activity. 

The remainder of this paper is organised as follows. In ^ we use 
neuroimaging data to construct a lattice-like network architecture em- 
bedded on the cortical surface of a rat brain, together with a sim- 
ple cellular automaton-like rule governing network activity. In ^ we 
present a comparative analysis of spreading dynamics over the cortical 
network, a uniform square lattice and an ensemble of 2D geometric ran- 
dom graphs. A summary of our main results, together with a discussion 
of future avenues of investigation, is provided in ^ 
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2. Spatial complex networks 



In this paper, we investigate the influence of network structure on 
the dynamics of the processes occurring upon them, with specific appli- 
cation to neural signalling. We achieve this by comparing the spreading 
dynamics of simulated neural activation within a rat cortex with those 
obtained on commonly-employed spatial networks. All of the networks 
examined herein are unweighted, undirected, and without loops. 

A plethora of tools and techniques for characterising complex net- 
works exists [31j; however, in the context of spatial networks not all of 
these measures remain relevant. An important feature of neural net- 
works is the propensity for nearby nodes to connect, and thus, such 
graphs tend to exhibit high levels of local clustering. In addition, ver- 
tex reachability (defined as the ability to travel between nodes i and 
j following connections within the network) impacts significantly on 
spreading dynamics. In view of this, systematic comparison of the 
network models considered here will be effected via the clustering co- 
efficient, which is given mathematically as 



where denotes the probability that any two neighbours of node i are 
connected and the order of the graph. Additionally, we employ the 
characteristic path-length, L, which is defined as the number of edges 
in the shortest path between two vertices, averaged over all pairs of 
vertices. We remark that if the network is not fully connected then the 
characteristic path-length diverg disconnected pair corresponds 

to an infinite path. 

2.1. Rat cortical network. Spatial coordinates defining the cortical 
surface of the rat were obtained from the Caret software package [32] 
and processed using the Caret Matlab toolbox. Figure [l|^a) shows the 
cortical surface of the left hemisphere of the rat brain, with typical 
neural activity spreading obtained via numerical simulation (see Q 
superimposed. 

Restricting to the left hemisphere, we construct a cortical rat net- 
work, with nodes positioned on the A^rat = 9623 available data points 
(see Figure [l](a)), by connecting nearest neighbours according to the 
following process: (i) a minimally-connected nearest-neighbour net- 
work is defined via the triangulation provided by the Caret software 
package; (ii) vertex pairs are connected if they lie within a Euclidean 
distance r of each other; (iii) connections are removed if the shortest 
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(a) (b) (c) 

Figure 1. A typical simulation of spreading dynamics 
in spatially constrained networks. Blue dots denote ac- 
tivated nodes (colour online), (a) A spatial network em- 
bedded on the left-hemisphere of the rat cortex; (b) a 2D 
square lattice graph; (c) a 2D periodic random geometric 
graph. 



path, calculated using the nearest-neighbour network defined in (i), 
between connected nodes exceeds a predefined number of steps, which 
is defined experimentally. While Euclidean distance is a reasonable 
estimate of physiological connection length in general [21], the third 
step is necessary to remove spurious edges which arise for large r, that 
are near in the ambient space, yet distant as measured on the cortical 
surface: shortest path length, defined via a simple mesh triangulation 
provides a simple method with which to measure this disparity. 



2.2. Standard network models. To discern the importance of net- 
work structure on activation spreading dynamics, we compare the spread- 
ing dynamics observed on the cortical network defined in ["lA with that 
observed on a uniform square lattice graph and on an ensemble of two- 
dimensional geometric random graphs [33], network structures typical 
of those employed in the literature [21] |22l |23j . 

To construct a random geometric graph we place uniformly and inde- 
pendently N nodes at random on the unit square, and form connections 
between pairs of nodes according to Euclidean distance. Similarly, a 
square lattice graph is obtained by forming distance-dependent con- 
nections within a A^-point uniform square lattice on the unit square. 
Figures [l]^b,c) illustrate these networks, and a typical simulation of 
neural activity. 

rt7^-^ we 



For comparability with the cortical network outlined in ^2.1 



choose A^ = 9604 (the closest square number to A'rat) and control the 
number of edges in the graph via f, the Euclidean distance scaled on 
the surface area of the rat cortex. The surface area is calculated via 
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the Caret Matlab toolbox to be ~ 221 (arbitrary units), which pro- 
vides a scahng: f = r/y/S. Additionally, we impose periodic boundary 
conditions. 

Figure [2] shows how the network measures C and L, with which we 
quantify differences between the three networks, vary with the connec- 
tivity distance, r. Ensemble measures for the random graph are calcu- 
lated from 10,000 realisations. These results indicate that the networks 
exhibit distinctly different features when restricted to short-range con- 
nections (small r), while highly-connected networks are comparable. It 
is noteworthy that the value of the mean clustering coefficient for the 
random graphs tends to the theoretical value 

(2) l-f(fc(i)'"°-'«*'^° 

obtained in p^, where r(a;) denotes the gamma function. 

We remark that for r < 0.2 we find f < 1/y/N, and therefore no 
connections exist in the square lattice graph (see Figure [2]^a)). At 
r = 0.2 the connectivity of the lattice graph corresponds to a four-point 
nearest-neighbour stencil (so that the connectivity matrix is analogous 
to a discrete Laplacian operator), and since nearest-neighbours of a 
vertex are not neighbours of each other, it follows that C = 0. In addi- 
tion, note that due to interplay between the regular discrete structure 
of the lattice and our rule based upon Euclidean distance with which 
to add new edges, variation of the clustering coefficient for the lattice 
graph is non-monotone with respect to r. The mean degree of the other 
networks is ^ 5 (data not shown). Lastly, we note that for r < 0.35, 
the random graphs are not fully-connected, hence the absence of data 
in Figure ^h). 

3. Spreading dynamics 



Propagation of activation within each neural network defined in §^2.1 



and 2.2 was governed by a basic spreading model [HI [IS], summarised 
as follows. 

Nodes i are restricted to exist in one of two states: active (xj = 1), 
or inactive (xj = 0). Starting from an initial activation state, simu- 
lation operated in discrete timesteps; from one timestep to the next, 
an inactive node became activated (or an active node remained in the 
active state) if it was connected to at least m active nodes. Initial con- 
ditions comprised a small region of activation (1% of the total nodes in 
the network) surrounding a node selected at random; the propagation 
of this activation through the network under our simple assumptions 
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(a) (b) 

Figure 2. Network measures. The clustering coefficient 
C and average path length L of each network plotted as a 
function of the connectivity distance, r. The dot-dashed 
line in (a) indicates the asymptote provided by Equation 
([2]). Data for the random graph is calculated from 10,000 
realisations; in the parameter range r G [0.35, 0.5] the 
maximum standard deviation of C and L is cr = 0.0016 
and 0.0301, respectively. 



on spreading dynamics provides a convenient and compelling method 
with which to highlight the differences imparted by the networks under 
consideration. We choose the mean fraction of activated nodes as our 
key metric with which to investigate the different networks. Ensemble 
measures of network dynamics on the random graphs and rat network 
were constructed from 10,000 simulations; behaviour in the uniform 
lattice is identical for all initial activation positions. 

The value of m places a lower bound on the connectivity of the net- 
work for which network activation can occur, and influences the speed 
of spreading of activation in the network (and, together with the value 
of r, the shape of the advancing activation front). Since we consider 
highly simplified dynamics in this study, omitting, for example, random 
inactivation or complex intra-node dynamics (the better to emphasise 
the importance of network structure on activation dynamics), the bal- 
ance between m and r determines completely the speed of activation 
of the network (indeed, for appropriate m and r, whole network exci- 
tation is inevitable) and, furthermore, affects all networks in the same 
manner. In all of the simulations that follow we chose m = 2 with- 
out loss of generality. For r < 0.35, the network structures are not 
comparable (as highlighted by Figure [2]^b)); for r > 0.5, the networks 
are highly-connected and activation spreads rapidly over the cortical 
surface. Our interest here is in the short-range connections between 
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cortical columns, rather than long-range white matter fibre tracts, say. 
We therefore restrict attention to connection distances r G [0.35,0.5], 
leading to mean degrees in the cortical network (rfrat) ^ [17.59,37.36]. 
Whilst actual estimates for mean degree in the rat cortex are not avail- 
able, such connectivity is typical of that employed in the literature (see, 
e.g. [9l [TOl [11] and references therein), particularly when one accounts 
for the large spread of observed degrees {e.g., for r = 0.35 the degree 
ranges from 8 to 48 whilst for r = 0.5 it lies between 20 and 87), and 
so serves to illustrate our methodology. For brevity, in the figures that 
follow, we illustrate the differences in spreading dynamics obtained in 
each network for the choices r = 0.35 and r = 0.5. 

Figure [T] shows a typical pattern of activation at an illustrative time- 
point in each of the three networks, which serves to highlight clearly 
how differences in the underlying network connectivity are made man- 
ifest in the spreading of activation through the network. We remark 
that patterns of excitation are 'well-defined' and uniform in the rat and 
lattice graphs, with all nodes within a region activated, while in the 
random graphs, lack of connectivity can lead to bottlenecks (here, we 
choose r = 0.25 to highlight these structural differences). 

Figure [3] shows the spread of activation through each network, as 
measured by the fraction of activated nodes (averaged over all in- 
stances). Figure [3]^a) indicates that for lower network connectivity 
(r = 0.35; the mean degree in each network was comparable: (rfrat) = 
17.59, ((irand) = 16.09 and diatt = 20) the activation speed in each 
of the three networks differs significantly: the lattice graph requires 
dramatically fewer timesteps to achieve entire network activation; the 
cortical network is the slowest of the three. In fact, we observe a 1.88- 
fold and 1.45-fold increase in activation time in the rat network when 
compared to the lattice and random graphs, respectively. For more 
highly-connected networks (r = 0.5), the cortical network remains the 
slowest to activate; however, the lattice and random graphs now dis- 
play similar activation rates, with the uniform lattice being marginally 
faster. The delay imparted by the cortical network is now 1.65- fold 
(lattice) and 1.46-fold (random graphs). 

Figures |4] and |5] highlight the infiuence of initial activation position on 
the spreading dynamics, which is quantified by the time to full-network 
activation, t* . For the values of connectivity distance, r, analysed here, 
full network spreading was observed for all networks independent of 
topology. Therefore, all simulations contributed to the calculation of 
t*. 

In Figures |4|(a)-(c) histo grams are presented, depicting the spread 
of t* observed in the simulations for each network; Figures |4Fd)-(f) 
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Figure 3. The evolution of the mean fraction of acti- 
vated nodes, p, in the rat network, random graphs and 
lattice for different connectivity distances: (a) r = 0.35, 
(b) r = 0.5. 



show the corresponding activation spread for each simulation, together 
with the ensemble mean. The delta-function obtained in Figure |4](b) 
and the corresponding results shown in Figure |4](e), reflect the fact 
that the network dynamics are identical for all realisations, due to the 
uniform structure throughout the lattice, and the periodic boundary 
conditions. In contrast. Figures |4]^a) and (d) indicate a very wide 
spread of activation times, and with no clear distribution, in the corti- 
cal network. In the random graphs, t* displays small variation on each 
realisation; a wider spread is exhibited for a less well-connected net- 
work (data not shown). As highlighted by Figure |4]^c), this distribution 
is well-characterised by a Gaussian with mean /x = 43.03 and variance 
a = 0.6901 (p < 0.01). The observed distributions of spreading dy- 
namics for the cortical and the random networks (Figures |4|^a) and 
[4](c) respectively) were found to differ significantly {p < 10^^ over the 
investigated parameter range) according to the Kolmogorov-Smirnov 
statistic. 

In Figure [5] heat maps are presented to highlight those initial acti- 
vation sites in the cortical network which provide significantly higher 
full-network activation speeds. These results indicate that the differ- 
ences in activation time shown in Figure |4] are due to the geometric 
structure of the rat brain and the presence of folded and smooth re- 
gions in the rat cortical surface. 

The results shown in Figure [5] are for the highly-connected network 
(r = 0.5); however, qualitatively similar results are obtained for the 
range of parameter values analysed here: although values of t* vary, in 
all cases, initial activation of the folded regions leads to (approximately) 
a 1.7- fold reduction in total of activation time (data not shown). 
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Figure 4. Simulation results indicating the variabil- 
ity in activation dynamics in each network. Panels (a), 
(d) contain results in the cortical network; (b), (e), the 
lattice; and (c), (f) the random graphs. Upper panels 
(a)-(c): Histograms showing the time to full activation 
t*. In (c), a comparison with a Gaussian distribution 
(yU = 43.03, cr = 0.6901) is also shown. Lower panels 
(d)-(f): The evolution of the mean fraction of activated 
nodes in each network (white circles), together with a 
confidence interval of width 2a. The connectivity dis- 
tance is chosen as r = 0.35. 



We have presented simulation results which highlight the differing 
rates of network activation, from an initial state comprising a small re- 
gion of activated nodes. Our model can therefore be thought of as 
broadly applicable to the initial stages of epileptic partial seizures, 
whereby spreading can initiate from a certain region, leading to in- 
creased activity patterns in larger parts of the brain (mechanisms that 
lead to the rise of initial activity, such as high frequency oscillations 
through gap junctions [33], or to seizure termination are not consid- 
ered). Additionally, the results that we have presented, indicating the 
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Figure 5. Heat maps indicating the dependence of ac- 
tivation dynamics on the initial activation site. The 
colours of each region indicate the time to full network 
activation, t*, associated with initial activation in that 
region (colour online). Panels (a) and (b) show rotated 
views; in each case r = 0.5. 



dependence on activation progression on initiation site are pertinent to 
observations in epilepsy, as not all brain regions have the same proba- 
bility of being starting points for seizures. We remark, however, that 
the concept of a specific focal point of seizure origin from which seizure 
activity spreads (due to local abnormal connectivity) has been criti- 
cised (see e.g. [36j), and alternative hypotheses presented j371|38]. We 
do not discuss this further as our focus here is on the global network 
structure imparted by the rat cortex embedding and its effect on acti- 
vation propagation. 

We remark that although our model has relevance to epileptic seizure 
spreading processes as described above, its form is intentionally simplis- 
tic: we have omitted a plethora of physiologically important network 
and signalling features (such as long-range connectivity, or complex 
node activation dynamics) as our study addresses a fundamental as- 
pect of the theoretical modelling of neural networks. Our central result 
is that cortical surface geometry affects profoundly signal propagation 
through the network, when compared to idealised networks with the 
same coarse statistics: the results contained in this section highlight 
clearly that significant differences in spreading dynamics are obtained 
in the network based on the cortical structure of the rat brain, com- 
pared to more standard network models. Since networks analogous to 
these standard models (defined in ^2.2) are employed frequently in the 
theoretical literature, in order to study the features of epileptic seizure 
initiation and dynamics, the simulation results highlighted by Figures 
|3| |4] and |5] indicate clearly that certain dynamics of relevance to seizure 
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initiation and progression (in particular, those highlighting the impor- 
tance of the site of initiation within the network) are not well-captured 
by such approaches. 

4. Discussion 

In this paper, we have investigated numerically the influence of the 
structure of a spatially-embedded complex network on the dynamics of 
the processes which occur upon it, with application to the development 
of improved theoretical models of the progression of epileptic seizures 
and spreading depression over the cortical surface. 

The key feature of this work is that wc employ neuroimaging data 
from the left hemisphere of a rat brain to define a neural network whose 
spatial embedding represents accurately the structure of the cortical 
surface. Typically, theoretical studies of neural network dynamics em- 
ploy uniform lattices or random graphs in ID or 2D: their focus being 
on in-depth analysis of various theoretical connectivity rules {e.g. small 
world or scale-free networks) or signalling processes on the network dy- 
namics. Here, we define connectivity within our cortical network by a 
simple criterion based on Euclidean distance (modified to account for 
the folded cortical structure), and employ a basic spreading rule to 
govern the propagation of activation. In this way, wc highlight clearly 
the influence of physiologically-relevant network structure on seizure 
dynamics, in isolation. 

Wc compared numerical simulation of the propagation of neural ac- 
tivity within three different network architectures: the cortical neu- 
ral network, a uniform square lattice graph and an ensemble of two- 
dimensional geometric random graphs, and studied these over a range 
of network connectivities. Note, that we limit our current study 
on comparisons of spreading on the brain surface and to short-range 
connections between nearby cortical columns; the role of long-distance 
flbre tracts between brain regions is not part of this work. In the 
parameter range chosen for the dynamic simulations which we investi- 
gate in detail, these networks were shown to be comparable in terms 
of their clustering coefficient and characteristic path length; in the 
case of networks with low connectivity {i.e. for values of r approaching 
the uniform lattice spacing) the networks display distinctly different 
characteristics. Despite their comparability, however, our simulations 
indicate dramatic differences in the propagation of activation through 
the various networks: in particular, we observe that the time taken to 
full activation in the cortical network is increased by a factor lying in 
the range 1.45-1.88. Especially striking is the importance of the site of 
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activity initiation: a dramatic spread in the number of time-steps re- 
quired to achieve full network activation is observed in the rat network, 
even in the case of a highly-connected network (r = 0.5, (c/rat) = 37.36). 
Small variation is observed in the random graph; the dynamics on the 
lattice graph is independent of initiation site. For less well-connected 
random graphs (r = 0.35, (errand) = 16.09), the variation is increased; 
however, the differences between all three networks remain significant 
(data not included). Moreover, our simulations highlight that, across 
the parameter range investigated here, different initial activation sites 
lead to approximately a 1.7- fold reduction in full- network activation 
time, due to variation in connectivity across these regions. 

We employ a highly simplified model with which to investigate the 
spreading of activation over our network. Various features of rele- 
vance to physiological neural networks, such as complex node dynamics, 
long-range connectivity, or stochasticity, are omitted from our formula- 
tion; indeed, our model draws no distinction between normal activation 
spreading, and that seen in epilepsy. However, such an approach al- 
lows us to provide a more powerful exposition of the importance of the 
network structure imparted by the cortical geometry, in isolation. We 
have shown that certain features of the activation dynamics display 
distinct differences, when compared to idealised models of the type 
commonly employed in the literature. Most strikingly, we highlight a 
variability in activation time for the cortical network, depending on 
initial activation site. This is pertinent to observations in epilepsy, as 
not all brain regions have the same probability of being starting points 
for seizures. Variability in cortical networks is well-studied in terms 
of structural and functional connectivity of brain regions, for exam- 
ple concerning the degree of nodes [391 SOI E] and its consequence 
on network robustness and performance |l2l |l3] . However, our result 
is remarkable in that variability is observed for curved brain surfaces 
even in the absence of white-matter fibre-tract connectivity. Our re- 
sults lead us to conclude that studies which do not take into account 
the spatial embedding of the cortex risk simplifying neural activation 
dynamics in a potentially significant way and are, therefore, unlikely 
to be able to represent accurately activation dynamics of relevance to 
epileptic seizures. We note, however, that in a physiological setting, 
similar disparity in network activation may be induced by including 
a range of other factors since the influences on the network dynamics 
within physiological neural networks are myriad; indeed, there is gen- 
eral agreement that no single factor can explain the varied phenomena 
associated with epileptic seizure dynamics. 
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We remark that epileptic seizure events are extremely rare in wild- 
type rodents, and characteristics of epileptic activity in animal models 
can differ from those observed in humans [35]. However, our initial 
study employing a network based upon the cortex of a healthy rat 
has highlighted the potential importance of geometric structure in ac- 
tivation dynamics; similar investigations in a network whose spatial 
embedding represents the convoluted human cortical surface therefore 
forms important ongoing work. In addition, future work will investigate 
how our predictions are altered under (i) a signal transmission model 
which represents more accurately the behaviour a neural unit and (ii) 
the introduction of long-range connections representing white-matter 
structural connectivity. 
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